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We present a comprehensive theoretical study of the static spin response in HgTe quantum wells, 
revealing distinctive behavior for the topologically nontrivial inverted structure. Most strikingly, the 
q = 0 (long-wave-length) spin susceptibility of the undoped topological-insulator system is constant 
and equal to the value found for the gapless Dirac-like structure, whereas the same quantity shows the 
typical decrease with increasing band gap in the normal-insulator regime. We discuss ramifications 
for the ordering of localized magnetic moments present in the quantum well, both in the insulating 
and electron-doped situations. The spin response of edge states is also considered, and we extract 
effective Lande ^-factors for the bulk and edge electrons. The variety of counter-intuitive spin- 
response properties revealed in our study arises from the system’s versatility in accessing situations 
where the charge-carrier dynamics can be governed by ordinary Schrodinger-type physics, mimics the 
behavior of chiral Dirac fermions, or reflects the material’s symmetry-protected topological order. 
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I. INTRODUCTION & MAIN RESULTS 

The two-dimensional (2D) topological insulator (TI) 
realized in inverted HgTe quantum wells exhibits un¬ 
usual electric-transport properties [1-4] that have at¬ 
tracted great interest [5]. TI behavior is also found in 
other 2D [6, 7] and bulk [5, 8] materials. The poten¬ 
tial for interesting interplays between a TI’s electronic 
and magnetic degrees of freedom, e.g., through hyper- 
fine interaction with the material’s nuclei [9], or local 
exchange interaction with magnetic dopants [10], has 
been pointed out recently [11, 12]. These efforts have 
opened up new perspectives and extended previous work 
devoted to understanding spin effects [13, 14] and mag¬ 
netism [15] in ordinary semiconductors. Our present the¬ 
oretical study of the spin response in HgTe quantum wells 
reveals unconventional spin-related properties that dis¬ 
tinguish this paradigmatic TI material from all other cur¬ 
rently known 2D electronic systems. We thus provide al¬ 
ternative means for the experimental identification of the 
topological regime and extend current knowledge about 
the fundamentals of spin-response behavior in solids. 

The spin susceptibility contains comprehensive infor¬ 
mation about the magnetic properties of a material. In 
the simplest case of a spin-rotationally invariant nonin¬ 
teracting electron gas, the spin susceptibility is propor¬ 
tional to the charge-response (Lindhard) function [16]. 
Noticable deviations from that situation occur, e.g., in 
systems with strong spin-orbit coupling such as 2D hole 
gases [17]. For metals or degenerately doped semiconduc¬ 
tors, the spin response of only the partially filled band 
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is typically considered. This approach misses intrinsic 
contributions to many-particle response functions aris¬ 
ing from virtual interband transitions that become im¬ 
portant in narrow-gap and, especially gapless, electron 
systems. Examples for the latter are the 2D Dirac-like 
electron states on the surface of a bulk TI whose mag¬ 
netic properties have been discussed in Refs. 18 and 19. 
The HgTe quantum wells considered here present an ideal 
testing ground for exploring the importance and proper¬ 
ties of intrinsic, or virtual-carrier, effects, as it is possi¬ 
ble to tune the band gap in such systems with a single 
structural parameter (i.e., the quantum-well width d). 
Furthermore, deviations from the conical 2D-Dirac dis¬ 
persion in the gapless case are well characterized within a 
continuum-model (BHZ) description [20] that also gives 
controlled access to the full range of, and interesting tran¬ 
sitions between, Schrodinger-physics-dominated, chiral- 
Dirac-fermion-like, and topologically nontrivial phases. 
In particular, definite (i.e., cut-off-independent) results 
for the intrinsic response at long wavelengths are ob¬ 
tained even in the limit of vanishing band gap; unlike 
in the case of the previously considered 2D-Dirac models 
used to describe the surface states of bulk TIs [18, 19]. 

Before going into greater detail in the remainder of 
this Article, we briefly highlight four major advances and 
central new insights gained from our work. 

(i) We find an analytical expression for the uniform 
static spin susceptibility of the intrinsic (undoped) sys¬ 
tem, 


xix*L(7;q = o) = 


clM 1 

167rlHl l + 4^M0(eM) 


( 1 ) 


where is proportional to the band gap and positive 
(negative) for the topologically trivial normal (topolog¬ 
ically nontrivial inverted) regime, 0 denotes the Heavi- 
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FIG. 1. (a) Mean-field critical temperature Tc(i) for virtual- 
carrier-mediated magnetism in a HgTe quantum well plot¬ 
ted as a function of the band-gap parameter ^m. Note the 
striking asymmetry between the topological (^m < 0) and 
normal (^m > 0) regimes. An explicit expression for To is 
given in the text; it contains the local exchange coupling be¬ 
tween magnetic-impurity orbitals and HgTe quantum-well ba¬ 
sis states, which is roughly constant across the transition [21]. 
(b) Critical temperature rc(d) for the electron-doped HgTe 
quantum well plotted as a function of and for various 
charge-carrier densities (corresponding to indicated values for 
the Fermi wave vector fcp; where fcp = qo corresponds to 
a charge density no = gQ/(27r) = 4.25 x 10^^ cm“^ for a 
typical structure [22]). We used = —0.7 for the BHZ- 
model electron-hole-asymmetry parameter in both plots, and 
7 = —2.22 as relevant for the (Hg,Cd,Mn)Te system [21]. 


side step function, S is a band-structure parameter from 
the BHZ model introduced further below, and the con¬ 
stants Cx^zil) depend both on the valence-band mixing 
of quantum-well basis states and on the parameter 7 
that characterises the relative coupling strength of the 
conduction-band and valence-band spin degree of free¬ 
dom to the physical quantity of interest. While the intrin¬ 
sic spin response is suppressed with increasing band gap 
in the normal regime (^m > 0 ), it becomes independent 
of gap size in the topological regime where < 0. This 
unexpected behavior is a special feature of the symmetry- 
protected topological phase of the bulk system associated 
with the existence of gapless edge states. 

(ii) A physical consequence of the unusual spin re¬ 
sponse in the intrinsic system is the asymmetric variation 
of the critical (Curie) temperature for virtual-carrier- 
mediated magnetic order in a HgTe quantum well that 


has been doped magnetically but not electronically. See 
Fig. 1(a). A similarly striking asymmetry arises when 
charge carriers are present in the 2D conduction band; 
which is illustrated in Fig. 1(b). Furthermore, in this sit¬ 
uation, a rather strong, and counter-intuitive, suppres¬ 
sion of the Curie temperature with the density of itiner¬ 
ant charge carriers is revealed. This tendency arises from 
the unconventional character of conduction-band states 
in the topological regime. See Section IIIB for details. 

(iii) We extract the effective g factors for 2D conduc¬ 
tion electrons in both the topological and normal regimes, 
which are directly accessible experimentally [4]. Their 
strong density and gap-parameter dependences reflect 
the band mixing and transition between regimes dom¬ 
inated by Schrbdinger-type and Dirac-like dynamics of 
charge carriers. Full details can be found in Section IV. 
These results are essential, e.g., to enable quantitative 
characterization of the transition between quantum spin- 
Hall and quantum anomalous Hall phases in magnetic 2D 
topological insulators [ 11 ]. 

(iv) We predict the effective g factors of the helical edge 
states that are present in the topological regime of the 
HgTe quantum well. For the out-of-plane magnetic-field 
direction, the g factor assumes a constant value that is 
determined by band mixing in the quantum-well eigen¬ 
states. In contrast, the in-plane g factor has a strong 
density dependence. See Section V. 

In the following, we discuss relevant details of our the¬ 
oretical analysis and present complete results for the spin 
response in both intrinsic and electron-doped systems. 


II. QUANTUM-WELL BANDSTRUCTURE 


The electronic properties of HgTe quantum wells are 
adequately captured by an effective four-band (BHZ) 
Hamiltonian [20] that acts in the low-energy subspace 
spanned by basis states |Ai-|-), \Hi+), \Ei—), |i?i—) and 
explicitly reads 


'n{k) 0 \ 

0 n*i-k)) 


( 2 ) 


with 'H(k) = h = (C — Dk'^, Akx, Aky, M — Bk'^) 

and = (1, cr). The parameters A, B, C, D, M are func¬ 
tions of the well width d, and their numerical values are 
typically obtained by a fitting procedure [23]. The pa¬ 
rameter M opens a band gap, where (in the convention 
B < 0) the system is in the topological (normal) regime 
when M < 0 (M > 0). The basis functions \Ei±.) are a 
superposition of conduction-electron and light-hole (LH) 
basis functions with a given spin projection. The asso¬ 
ciated band |Ti±) has a much lower band-edge energy 
and can therefore be omitted. On the other hand, the 
heavy-hole (HH) states \Hi±) also belong to the set of 
low-energy excitations. In the following we set C = 0, 
and employ a dimensionless description of Eq. (2) that 
is obtained by defining an energy scale Eq = AqQ and a 






3 


scale for the wave vector = A/\B\ [24, 25]. (In a typi¬ 
cal HgTe quantum-well structure [22], Eq = 0.189 eV and 
go = 0.517 nm“^.) By making use of the axial symmetry 
of the BHZ Hamiltonian (2) we rotate to a real basis to 
obtain 7^'(k) = = Eoh’^a>^, with = 

diag(e**<^'‘/2,e-**‘^>'/2) ^nd E = + k^), 

where 0k is the polar angle of the 2D wave vector k. We 
have defined the dimensionless parameters = M/Eq 
and = D/\B\, which have typical values [20, 22, 23] 
I^mMCd] ^ 0.5. The eigenvectors in the complex and 
real basis are related via ®iai where a = ± 

distinguishes the conduction and valence bands, which 
are doubly degenerate in the quantum number s = ± for 
spin projection along the growth direction and have the 
dispersions 

4a ^ 4a = Eo[-~k^^j, + a^{U + ~BY+k^] . (3) 


III. EFFECTIVE SPIN SUSCEPTIBILITY OF 
THE 2D SYSTEM: INTRINSIC & DOPED CASES 


The spin susceptibility is used to characterize a sys¬ 
tem’s response to spin-related external stimuli within the 
framework of linear-response theory [IG]. In the static 
limit, and for quantum-well-confined electrons, it can be 
written as [17] 


Xy(R,z;R',z') = lim 




dt e-’'* 


x([5,(R,z;t), 5,(R',z';0)]) 


, (4) 


therefore need to express the spin susceptibility within 
the underlying six-band Kane model. 

The coupling between some physical stimulus repre¬ 
sented by a field and the Tg and Tg-band intrinsic 
angular-momentum degrees of freedom a and J is most 
generally described by a term 

® ^4x4 + 02x2 ® Xj (6) 

in the Kane-model Hamiltonian. See, e.g.. Table C.5 in 
Ref. 26. Within this approach, the coefficients br are 
intra-band coupling constants with appropriately renor¬ 
malized values to take account of all field-induced band¬ 
coupling effects in the bulk material. To be able to dis¬ 
cuss a wide range of spin-related phenomena, we define 
an effective (pseudo-)spin operator 




(7) 


such that = 6 rg J2i Ei Si{brgjbrs)- The actual value 
of the parameter 7 depends on the physical quantity or 
situation of interest [27]. The system’s response is then 
fully captured by the effective spin-susceptibility tensor 





'^zj'(k,k+q,a,^) 


ir^z.z') 




(s’) 


ihr] 


e21-e. 


{s') 

k-l-q^ 


( 8 a) 


where np denotes the Fermi function, and 


with R = (a;, y) and z being coordinates in the 2D plane 
and perpendicular to it, respectively, and ^^^R, z) = 
4'f(R, z) Sj 4'(R, z) denoting the electron spin density 
measured in units of ft. For a homogeneous 2D electron 
system, the spin susceptibility is most straightforwardly 
obtained in terms of the spatially Fourier-transformed 
quantity Xij(q; z, z') via 

X.,(R,z;R',z')= I (5) 


^zj'(k,k+q,a,^) 


{r,z,z')=ki^^{z) 


Sii-j) 4+4 (^) 


4 + 4 ^ 


it 


(z') ' • 5j(7)4a(^') 


( 8 b) 


are matrix elements of the effective spin operator given 
in Eq. (7). In the spirit of subband k • p theory, the 
six-dimensional spinor wave functions (z) can be ex¬ 
pressed in terms of the BHZ-model basis state spinors 
44 ^) as 


The dependence of Xy(q;z,z') on the coordinates 
z and z' encodes the spatial profile of the quantum- 
well bound states. Within the BHZ framework, the z- 
dependent part of electron wave functions is contained 
in the four basis functions The latter are 

spinors whose explicit form has been derived [ 20 ] within 
the six-band Kane-model description [26] for the charge- 
carrier dynamics that includes the bands with Fg and 
Fg symmetry closest to the bulk-material’s fundamen¬ 
tal gap. As is generally the case in multi-band sys¬ 
tems, the spin response of electrons in a HgTe quantum 
well is strongly influenced by both the in-plane dynam¬ 
ics described by the BHZ Hamiltonian and the nontriv¬ 
ial spinor structure of the BHZ-model basis states. We 


44^) = E(4i0.4S*^4W: (9) 

i=l “ 

fs) 

where the coefficients ^ are the components of the 
corresponding eigenvectors of the BHZ Hamiltonian (see 
Sec. H). The explicit form of the basis states was derived 
in Ref. 20 by solving a confined-particle problem in the 
HgTe/CdTe hybrid system. For instance, 044)^ = 
(/i(z), 0 , 0 ,/ 4 (z), 0 , 0 ), with the detailed form of the en¬ 
velope function components fi{z) provided in the supple¬ 
mental information of Ref. 20. 

In the following, we consider the growth-direction- 
averaged spin susceptibility of charge carriers in 
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the HgTe quantum well, given by Xy(7;q) = 

f dzf dz' [17], with z, z') calcu¬ 

lated using the Kane-model-based BHZ approach as de¬ 
scribed above. Note that, by replacing the spin matrices 
in Eq. (8b) by the unit matrix and averaging over z, z\ we 
obtain the charge-response function studied in Refs. 24 
and 25. From the axial symmetry of our Kane-model 
description, and the fact that the eigenstates have def¬ 
inite spin projection in the growth direction, it follows 
that the in-plane spin susceptibilities are the same, i.e., 
Xxxil: q) = Xyyp', q), and only depend on the magnitude 
q = |q|. This is an important difference to other semi¬ 
conducting systems where HH-LH mixing occurs [17]. In 
the present situation, HH-LH mixing arises only from 
terms linear in k (giving rise to Dirac-like excitations) 
which is a consequence of the envelope function compo¬ 
nents [20] behaving differently (even or odd) under the 
parity transformation z -O- — z. 


A. Spin response of the intrinsic system 

In the insulating limit, the conduction band is empty 
and the Fermi level lies in the band gap, i.e., \fi\ < 
T^oI^m]- In this situation, the spin susceptibility origi¬ 
nates from virtual interband transitions across the band 
gap. In the following we will focus on the zero temper¬ 
ature limit. The analytical result for the intrinsic spin 
susceptibility obtained in the limit of zero momentum 
transfer g —>■ 0 is given in Eq. (1), where 

Cx(7) = 1 + (27 - 1)Clh , (lOa) 

(^ 2 ( 7 ) = 1 - 37 -I -(7 - I)Clh , (lOb) 

and Clh = / dz\f 4 {z)\‘^ (~ 0.4) is the amount of LH 
admixture in the basis states lFi±). (Details about the 
calculation of the intrinsic spin susceptibility are given in 
Appendix A. Results for g 0 can be calculated numer¬ 
ically. For completeness, some of these are also shown in 



FIG. 2. Prefactors and Cy determining the magnitude of 
intrinsic in-plane and out-of-plane spin response, respectively 
[see Eq. (1)], plotted as a function of the parameter 7 from 
the definition of the effective Kane-model spin operator (7). 


Appendix A.) Equation (1) exhibits anomalous behav¬ 
ior in the inverted region (^m < 0 ) in that the uniform 
spin response is independent of the gap size and pinned 
to the value for the gapless case, even though it arises 
from virtual inter-band transitions. In the normal region 
(^M > 0 ), the expected decrease of the response functions 
with increasing band gap is found. The spin suscepti¬ 
bility ( 1 ) is strongly anisotropic, generally exhibiting a 
dominant out-of-plane response except for a very small 
range of the parameter 7 . See Fig. 2. 

A non-vanishing ( 7 J 0) seems to imply the coun¬ 
terintuitive phenomenon that an applied magnetic field 
could generate a magnetisation of the HgTe quantum- 
well system in the intrinsic limit where no charge carri¬ 
ers are present. Our more detailed analysis shows, how¬ 
ever, that this is not the case. Direct calculation of the 
derivative of the free energy with respect to magnetic- 
field strength in a model based on the BHZ Hamiltonian 
(2) augmented by a Zeeman term reveals that the total 
magnetisation is strictly zero [28]. This is due to the 
fact that intrinsic-spin and orbital contributions to the 
magnetization cancel, as expected in a spin-orbit-coupled 
system [29]. This conclusion is further underpinned by 
the observation that Xj 7 *^^( 7 ; 0 ) vanishes in the limit of 
zero HH-LH mixing [30]. 

As it is possible to engineer and study effective ex¬ 
change interactions between impurity atoms [31, 32], it 
is tempting to consider such interactions between two lo¬ 
calized spins in a HgTe quantum well. Quite generally, 
the RKKY Hamiltonian is given by [33] 

^r':®=G2E ^r7^^Xb(r,r'). (H) 

Here Xy (r, Fourier transform of ( 8 a), is the spin 

susceptibility in real space, I® denotes the ith Carte¬ 
sian component of an impurity spin located at position 
r = (R, z), and G is the local exchange-coupling constant 
between the spin degree of freedom carried by band elec¬ 
trons and localized (e.g., impurity) spins. (The difference 
in exchange-coupling strengths for conduction-band and 
valence-band states is accounted for by the appropriate 
value of 7 .) In Fig. 3(a), we plot the growth-direction- 
averaged effective out-of-plane real-space spin suscepti¬ 
bility Xzz(77) as a function of the distance R for various 
values of (and = —0.7 corresponding to realistic 
situations [23]). 

Figure 3 shows that the effective exchange interaction 
mediated by the intrinsic spin response of the HgTe quan¬ 
tum well for localized spins is of ferromagnetic (FM) type 
if Ago ^ I, and that there is a cross-over to antiferro¬ 
magnetic (AFM) coupling for Ago > I that sensitively 
depends on ^m- In particular, we find that this cross¬ 
over is shifted significantly to smaller impurity-spin sep¬ 
arations in the case of the inverted regime as compared 
to the normal regime. Also, the magnitude of exchange 
interaction in the ferromagnetic regime for fixed distance 
R markedly increases when enlarging the gap in the in¬ 
verted system. In contrast, the exchange interactions in 




5 



FIG. 3. (a) Intrinsic real-space spin susceptibility Xzzi^) 

plotted as a function of Rqo for several values of the gap 
parameter ^m, with fixed = —0.7 and 7 = —2.22. (b) Pa¬ 
rameter domains of preferential alignments between localized 
Mn spins having a distance R, mediated by the intrinsic spin 
response of an HgTe quantum-well system. The blue (white) 
domain favours a ferromagnetic (anti-ferromagnetic) out-of¬ 
plane spin alignment. In the yellow region the two localized 
spins tend to align ferromagnetically in-plane. 


the normal regime have a very small magnitude. To fur¬ 
ther illustrate the parametric dependencies of the carrier- 
mediated exchange interaction, Fig. 3(b) shows the var¬ 
ious regions of preferential spin alignments as a func¬ 
tion of < 0. Interestingly, inbetween the two out-of¬ 
plane localized-spin alignments (FMj_ and AFM^), we 
also find sizeable regions for the distance where it is en¬ 
ergetically favorable for impurity spins to align ferromag¬ 
netically in-plane (FM||). Moreover, for larger distances, 
we find the usual Bloembergen-Rowland [34] behavior for 
the intrinsic spin susceptibilities which is characterized 
by an exponential decay factor that depends on the band 
gap, approximately given by ^ exp {—R/Rq)/ where 
^0 = (ICmI^o)”^ is the Compton wave length of the band 
electrons. More details are given in Appendix B. In the 
limit of zero gap (^m —>■ 0 ), the R~^ decay found previ¬ 
ously in various Dirac systems [18, 19, 35] is reproduced. 

Having considered the case of two localized impurity 
spins with exchange interactions mediated by virtual ex¬ 
citations, we now focus on a HgTe quantum well that is 
doped with a large number of homogeneously distributed 

I 


magnetic impurities (eg, Mn ions), forming effectively a 
Hgi_ 3 ,Mna;Te alloy [36]. The average distance between 
the spins is i?NN ~ (3ag/167ra;)^/^ [37], with gq the HgTe 
lattice constant and x the concentration of magnetic ions. 
We calculate the Curie temperature in the mean field 
limit, assuming R^jsi/Rq 1 which is justified when 
X :$> 2 X where we have used qo ^ 0.5 nm“^. 

Since j.jM| ^ 0.3 [23], this condition is practically al¬ 
ways fulfilled. The Curie temperature for Ising-type fer¬ 
romagnetic order with magnetization perpendicular to 
the quantum well growth direction is given by 


- , dc/dj^u) 

®0) 0 1 + 4^m0(Cm) 


( 12 ) 


where Tg = ^ denotes the impurity- 

spin magnitude, nj is the 3D density of magnetic im¬ 
purities, and dc ~ 6.3 mn is the critical well width. 
For obtaining Eq. (12) we have used the approxima¬ 
tion /-f/2 (7; 0, z,z') Ri (o'’ 


Fig. 1(a), we show the Curie temperature as a func¬ 
tion of 5 m, where we set A = 0.375 eV nm and B = 
— 1.120 eV nm^ [23] as their dependence on d is much 
weaker compared to M. We see that the behavior of the 
Curie temperature in its dependence on 5 m, or equiva¬ 
lently d, provides a clear means to distinguish between 
topological and normal regions. 

Knowing the spin susceptibility as a function of the 
wave vector allows us to go beyond the mean field 
limit [13, 14] and discuss the stability of the mean-field 
ground state with respect to thermally excited spin waves 
(magnons). In our present case of interest, the magnon 
dispersion for g/qg -C 1 is given by ujq = wq + C 2 (?^, with a 
1 


coefficient C 2 = 


dq^ 


q^O 


. Using Figure 6 in Ap¬ 


pendix A, which shows results for the relevant set of pa¬ 
rameters, we find C 2 < 0. Thus the criterion [14, 17, 33] 
C 2 > 0 needed to guarantee stability of the mean-field 
ground state is generally violated. It may be possible that 
electron-electron interactions and/or spin-orbit coupling 
help to stabilize ferromagnetic order, as has been shown 
to be the case for an ordinary 2D electron gas [13, 14, 38], 
but this question is beyond the scope of our present work. 


B. Spin response of the electron-doped system 

We consider now situations where the conduction band 
is occupied (/i > Uo15m|). In the g —^ 0 limit, an analyt¬ 
ical expression for the extrinsic contribution to the spin 
susceptibility (arising from filled conduction-band states) 
is found as 
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v(dop) 

\xx 


(7;0) 


TfePi —- 2 /(% + {m) - 2 / j +-TiTTlfS- 


(13a) 


Cl{^) r 27rjV(0) / 67P.(7) p ^f{Vl{lM+^M + fY + <^l^~4} \ 

IQtt\B\[ P \ C2(7) f ci{p[kl + {kl+^M + f?Y ) 


fcp(l + 2 ^m) + 2^^ — 2/|^m|) 

/(I + 4^m) 


(13b) 


where the first (second) term between the square brackets 
in the expression for Xj 7 °^^( 7 ; 0 ) is an intraband (inter¬ 
band) contribution. In Eqs. (13), we used the abbrevia¬ 


tions ( 7 ) = (1-7)Clh-1 and / = \J (fc| -|- ^m)^ + 
iV(0) is related to the density of states at the Fermi en¬ 
ergy as -Xo(0) = N{Q) = {2/\B\) iV(0) (where Xoio) is 
the static charge-response function), and reads -/V(O) = 


1 

277 


l~l~ 2 /cp + 2 ^m 


— 2 ^d 


For the general case with 


"x/ (il‘F+^M)^-t-fcp 

9 > 0 , Xjjin'iO) can be calculated numerically and its 
line shape is presented in Appendix C for various ex¬ 
amples. Due to the sharpness of the Fermi surface, we 
find for the corresponding real-space spin susceptibilities 
the expected R~^ oscillatory decay expected from a 2D 
Fermi liquid [16]. 


Next we calculate the Curie temperature using the 
mean field approach under the premise ^ 1 , 

which leads to the condition x ^ 2 x 10“^ k-p- This 
condition is generally fulfilled since kp < 0.1 nm“^ is 
the reliable range of the effective model [23]. The Curie 
temperature is given by 


Tc(d) = To 


dc 167r|i3| 

^ 'cM 


lxz^(7;0)| 


(14) 


where X^^( 7 ; 0 ) = X&*^(7;0) + xS°‘’^(7;0)- Note that 
also in the doped case Ising type magnetism prevails. In 
Fig. 1(b), we show the Curie temperature as a function 
of /m for various values of the Fermi wave vector. For 
low doping, we see a strong dependence of the Curie tem¬ 
perature on in the topological regime, which becomes 
very weak in the normal regime. For large doping, on 
the other hand, the Curie temperature is not very sen¬ 
sitive to in its whole range. Furthermore, the Curie 
temperature is generally suppressed with increased dop¬ 
ing, but this trend is much stronger in the topological 
regime. For the doped case, we find the magnon disper¬ 
sion ujg = coo + C 2 q^ where C 2 = ^ ^ ^ 

^ 9=0 

always, in contrast to the intrinsic case. See Fig. 8 in 
Appendix C. Thus a necessary condition for the stability 
of the ferromagnetic groundstate is fulfilled [14, 17, 33]. 


IV. CONFINEMENT DEPENDENCE OF THE 
2D-ELECTRON EFFECTIVE g FACTOR 


The paramagnetic response of charge carriers in 
quantum-confined structures is usually interpreted in 
terms of the Pauli spin susceptibility and quantified by 
an effective single-particle g factor [39-41]. However, any 
actually measured spin-related quantities correspond al¬ 
most always to averaged collective responses of the, e.g., 
quasi-2D, electron system that can be crucially affected 
by nontrivial spin-related phenomena [17]. Here we con¬ 
sider the paramagnetic response of conduction-band elec¬ 
trons in HgTe quantum wells and show how their param¬ 
agnetic response is changed as a function of the band-gap 
parameter that drives the transition between the topo¬ 
logical and normal regimes. 

To define an effective g-factor for our system of a HgTe 
quantum well, we introduce the bulk-material Zeeman 
term Jfs = 9 * Sj{—2K /where /ie is the 

Bohr magneton and g* (k) the Tg-band (Tg-band) g fac¬ 
tor [26]. Linear-response theory enables us to deter¬ 
mine the paramagnetic response to the magnetic field, 
which is given by Xpj = ( 5 *Mb)^ Xjj°‘’^(- 2 k/ 5 *;'?)l<?=o, 

where ^^°^\—2Klgj,\q) are the spin susceptibilities of 
the electron doped system for the in-plane and out-of- 
plane response involving Sj[—2K/gj^), see Fq. (13) for 
7 = —2K/g*. We compare this with the Pauli suscepti¬ 
bility given by Xp^ = {gjmf Xo(0)/4, with Xo(0) being 
(up to a minus sign) the density of states which is the 
zero-g limit of the Lindhard function and gj is the Lande 
g-factor for the two directions. Thus, we can extract col¬ 
lective g factors for the charge carriers as 


9j — 9* 


\ 




Xo{q) 


q=0 


(15) 


Our approach to determine g factors via the spin suscep¬ 
tibility complements previous work [42] where an effec¬ 
tive Zeeman term was derived for the BHZ Hamiltonian. 

For a narrow-gap semiconductor, the values for g, and 
K are dominated by band-coupling contributions. Using 
generic expressions resulting from Kane-model descrip¬ 
tions [26], we find 


2 k ^ + Ca 

9* 2^a 


(16) 


with = Aq \B\I in terms of the spin-orbit-splitting 
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gap Aq between the Fg and Fr band edges in the 8x8 
Kane model. As ^ |Cm| for our situations of interest, 
we can set 2k/ g^: —>■ —0.5 in the following. Hence, in 
Fig. 4, we show the effective ^-factors for in-plane and 
out-of-plane responses as a function of for various 
levels of electron doping obtained for 2 k,/ = —0.5. 

For very low densities, we see a vanishing in-plane re¬ 
sponse and a maximal g^ = 6 k HH-like out-of-plane re¬ 
sponse in the topological region < 0. This behavior 
arises because the conduction-band character is domi¬ 
nated by the HH basis states, which experience a “frozen” 
spin orientation perpendicular to the quantum well due 
to the confinement-induced HH-LH energy splitting [43] . 
This can be easily verified from (13) by taking the limit 
k-p 0. Departures from these results occur for larger 
doping levels (larger kp) due to increased HH-LH mixing. 
In the normal region (^m > 0), we encounter the situa¬ 
tion that the conduction band is dominantly composed 
of the electron basis states, rendering the g factor to be 
sizable at any doping. The fact that the in-plane spin 
response is notably larger than the out-of-plane response 
for ■fM > 0 and low doping is due to the LH admixture 
in the conduction-band states. As the carrier density in¬ 
creases, the concomitantly increased HH-LH mixing re¬ 
sults in significant modifications. More detailed explo¬ 
ration of parametric dependences exhibited by the col- 


_l ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 

(a) 

^ kF=0.3qQ —■ 

._ 



l<F=0.5qo _ 

/ 

/ 

^F=0.01(7o ^ 

.... 1 .... I .... i~ 


-0.3 -0.2 -0.1 0 0.1 0.2 0.3 



-0.3 -0.2 -0.1 0 0.1 0.2 0.3 


FIG. 4. Effective p-factors associated with (a) the in-plane 
and (b) the out-of-plane response of an electron-doped HgTe 
quantum well, plotted as a function of the gap parameter 
for various levels of doping and with = —0.7, 2k/ gt = 
—0.5. (Note that fcp = go corresponds to a charge density 
no = go/(27r) = 4.25 x 10^^ cm“^ for a typical structure [22].) 


lective g factors obtained here will be useful to augment 
previous perturbative estimates [ 12] and aid the interpre¬ 
tation of recent measurements [4]. 


V. SPIN RESPONSE OF QUASI-ID HELICAL 
EDGE STATES IN THE TOPOLOGICAL REGIME 

We have also investigated the edge-state contributions 
to the spin susceptibility in the topological region, finding 
them to be negligible compared to the bulk contributions 
in almost all situations. The only exception occurs for the 
in-plane response function Xxxin ?) the purely intrin¬ 
sic situation where the Fermi level lies in the minigap of 
the edge-state dispersions that opens up in a finite-size 
sample [5, 44]. More details and full results pertaining 
to edge states are given in Appendices D and E. 

Upto very small finite-size corrections, the paramag¬ 
netic response of helical edge states is captured by using 

(e) 

the effective ^-factors g//z given by 



(17b) 

The results shown in Eqs. (17) were obtained using 
Eq. (15) together with Eqs. (E7a)-(E7b). The prediction 
of the edge-state g factor is a major result of our work. 
In Fig. 5, we plot the in-plane g-factor as a function of 
kp where we used the natural cut-off scale A = tt/oq. In¬ 
terestingly, the in-plane g-factor of edge states decreases 
monotonically with increasing doping level. Such a be¬ 
havior is in stark contrast to the bulk case shown in 
Fig. 4(a) and to other 2D systems [17, 45] where the 
in-plane response increases with increased doping. This 
behavior reflects the fact that the spin-quantisation axis 



FIG. 5. Effective g-factor characterising the response of heli¬ 
cal edge states to an in-plane magnetic field. 





















of the edge states is perpendicular to the 2D plane, i.e., 
parallel to the 2 direction, and their helical nature. 


VI. CONCLUSIONS 

We present a detailed theoretical study of the spin 
response in HgTe quantum wells where virtual-carrier- 
related processes are particularly relevant and exactly 
tractable within the effective BHZ-model description. 
Anomalous properties of the spin susceptibility and 
carrier-mediated magnetism are found in the inverted 
regime, extending our current understanding of spin- 
related properties in topological materials. 

Most strikingly, the uniform static spin susceptibility 
in the intrinsic limit is constant, independent of the band 
gap, in the topological regime [see Eq. (1)]. This results 
in a distinct asymmetry between normal and inverted 
systems illustrated, e.g., by the dependence of the critical 
temperature for virtual-carrier-mediated magnetic order 
in a system that is doped with magnetic ions. Impor¬ 
tant differences between the two regimes are also exhib¬ 
ited in the situation where the conduction band becomes 
filled (see Fig. 1). The stability of mean-field ferromag¬ 
netic ground states with respect to thermal excitations of 
magnons has been analysed. We find that the magnetic 
order is stable in the situation with finite doping. 

In the topological regime, quasi-ID helical edge states 
exist. We have investigated their spin-response proper¬ 
ties, finding that their contribution to the spin suscepti¬ 
bility is thermodynamically suppressed compared to that 
arising from 2D quantum-well states, except in the very 
special - and probably physically hard-to-realise - sit¬ 
uation when the chemical potential is in the mini-gap 
opened by the hybridisation of states from opposite edges 
in a finite sample. 

We have used our results obtained for the spin suscep¬ 
tibility to define effective collective g factors for states 
from the quasi-2D quantum-well subbands and also for 
the quasi-ID helical edge states. The different character 
of quasi-2D-subband states in the normal and topological 
(inverted) regimes is reflected in the values for the effec¬ 
tive g factor. Their strong dependence on charge-carrier 
density reveals the importance of interband mixing. The 
behavior of the g factors found for the edge states reflects 
their helical nature and spin-quantization property. 

Our results are directly relevant for current and poten¬ 
tial future experimental investigations of the spin-related 
properties of 2D topological insulators, in particular 
those realised in HgTe/HgCdTe and InAs/GaSb quan¬ 
tum wells. For example, recent observation of Josephson- 
junction interference patterns in an S-HgTe/HgCdTe-S 
hybrid system has enabled extraction of g factors for the 
quantum-well charge carriers [4]. It would be interesting 
to use similar techniques [46, 47] to measure the edge- 
state g factors and compare with our predictions. Fur¬ 
thermore, our results for the spin response in both the 
intrinsic and doped regimes are informative for the design 


of, and interpretation of measured quantities for, dilute 
magnetic phases in these systems [10, 15]. 

It would be interesting to extend our formalism to 
study the spin response in 3D topological-insulator ma¬ 
terials [8]. In particular, as was the case in the 2D 
quantum-well-based TIs considered in this work, the in¬ 
terplay of charge-carrier dynamics and the spinor char¬ 
acter of extended bulk states could be a source of rich 
variety in spin-related properties also in 3D, and the con¬ 
tributions of the conducting surfaces are currently not 
understood. 
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Appendix A: Intrinsic contribution to the spin 
susceptibility 


The intrinsic contributions to the diagonal entries of 
the spin susceptibility (the only ones that are non-zero 
for the BHZ model) are calculated by 

-(int)/ N ^ f (Pk ^j(’k,k+q, + 

^ £^k^ - ^k+q- + ik.r]S 

5 = ±1 ^ 

(Al) 


where in the zero-temperature limit the valence band 
is fully occupied, i.e., nF(F^'!l) = 1. Here and 

in the following we consider the growth-direction- 
averaged case: Xy ( 7 ;q) = J dz J dz' Xy ( 7 ;q;z,z') ^ 




W 

Wi(k,k+q,a,/3) 


( 7 ) = jdz dz' In 

Fig. 6, we plot Uxxil'^d)] in panel (a) [(b)] 

as a function of q for various values of ^m- The inset 

< 0 for all 


Fig. 6(a) illustrates that ^ 


9=0 


of the values of considered, which implies that the 
mean-field ferromagnetic order for out-of-plane aligned 
magnetic impurity spins will generally be destroyed by 
spin-wave (magnon) excitations [14]. 
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FIG. 6. Intrinsic spin susceptibilities (a) Xzzil’^) (b) 
Xxxi'y':^) function of q/qo for various values of = 

(—0.3,—0.2,—0.1, 0,0.1, 0.2, 0.3) (from top to bottom) with 
= —0.7, Clh = 0.4 and 7 = — 2 . 22 . The inset in (a) 
displays the normalized spin susceptibility XzziTi <l)/Xzz{'y'i 0 ) 
in the small q/qo region for = (0.1, 0.2, 0.3), which shows 
more clearly that Xzz (7! <?) > Xzz (7! 0) iu the small-g limit 
also in the normal regime. 


Appendix B: Bloembergen-Rowland behavior of the 
local intrinsic spin susceptibility 


Bloembergen and Rowland [34] found that the local 
RKKY interaction of gapped systems becomes short- 
ranged, i.e., is exponentially suppressed by the band-gap. 
The functional dependence of the local spin susceptibil¬ 
ity on the distance for the gapped Dirac system at hand 
can thus be modelled by 


Xjj{R) ~ 



(Bl) 


where = |^m| 9 o is the inverse of the Compton wave¬ 
length of the system and 0 ( 1 ) is a numerical coef- 
hcient that can depend on the distance itself. To illus¬ 
trate this behavior for the HgTe quantum well system, we 
plot Xzz{R) ill Big- 7 as a function of Rqq for = 0.1 
[^M = — 0 . 1 ] in panel (a) [(b)] together with both the line 
shape expected for 2D massless-Dirac particles and the 
Bloembergen-Rowland result. 




FIG. 7. Local intrinsic spin susceptibilities Xzz{R) (solid or¬ 
ange curve) as a function of Rqo for (a) =0.1 and (b) 

= —0.1 with = —0.7, Clh = 0.4 and 7 = — 2 . 22 . The 
dashed green curve shows the modelled behavior in Eq. (Bl) 
with coefficients c = 1 for (a) and c = 1.2 for (b). For compar¬ 
ison we also show the R~^ decay for a gapless Dirac system 
represented by the dot-dashed blue curve. The sharp drop of 
the orange curve in (a) is due to a sign change of XzziR) 
about Rqo ~ 8 . 


Appendix C: Electron doped contribution to the 
spin susceptibility 


For the case where the Fermi energy is above the con¬ 
duction energy band edge, i.e., fx > \M\, the spin sus¬ 
ceptibility receives contributions due to electron doping 
given by 


X 


(dop) 


( 7 ; q) = XI 


(Pk 






Bk+ - Bk+q+ 


ihrjS 


^!;^k2+q,+,-)(7) ' 


(Cl) 

where for zero temperature np(7?^+) “ ® ~ I^^I) ; with 

fcp being the Fermi wave vector associated with the con¬ 
duction band. The complete spin susceptibility in the 
doped case is therefore given by 

Xjj ( 7 ; q) = Xj”*^ ( 7 ; q) + ( 7 ; q) • (C2) 

We show the line shape otXzziXT^l) and g) in Fig. 8 

for = ±0.2 and various choices of doping. The insets 
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FIG. 8 . Normalized spin susceptibilities Xz 2 ( 7 i l) Xa:a;(7) l) (b,d) as a function of q/{2kp) for varions levels of doping. 

(a,b) are the results for = — 0.2 and (c,d) for = 0 . 2 , with = —0.7, Clh = 0.4 and 7 = — 2 . 22 . The insets of (a,c) show 
the qnadratic q = q/qo dependence close to g = 0 . 


of Fig. 8(a,c) show the small q dependence of Xz2(7;9); 
indicating the stability of the out-of-plane Ising-type fer¬ 
romagnetism with respect to Magnon excitations. 


can be written as 


k,Q;,s 


gZk-R 

/z7 


V'ka(e) C 


(s) 

kcn 


(D2) 


Appendix D: Combining bulk and edge 
contributions to the spin susceptibility 


In order to compare bulk and edge contributions to the 
spin susceptibility we begin by considering the real-space 
spin-response function 


7 

Xy (r,r') = -- J dt e {[S,{v,t), ^^(r^O)]) , 

(Dl) 

where ^^(r) = dff(r)^ilE'(r) are spin-density operators. 
In our system of interest, electrons shall be confined to 
move freely in d < 3 dimensions, with system size L in 
all of these free directions. The position vector shall be 
split up into a part R comprising the coordinate direc¬ 
tions in which the motion is free and a part q in whose 
coordinates motion is confined; r = (R, q). The second- 
quantised electron operator in real-space representation 


with normalized spinor bound-state wave functions 
Cnk(^)- A straightforward calculation yields 


X^3\ 


= :^ E Q, q') , (D3a) 


with the q-dependent spin susceptibility of the d- 
dimensional (dD) system given by 


fdD)/ 


Xy ^(q;e,eO= E E ^S(k.k+q.a,/3) 

0',^,s,s' k 




^k+ql3 


■ ihr] 


.(D3b) 


We are now interested in the homogeneous part of the 
spin response defined as 


T — 


d^r 




r x^3{r,r')^L^tr{^ = 0) , 

(D4a) 
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where 

= J d^~‘^g J £>,£>') ■ 

(D4b) 

In the situation where both 2D bulk and ID edge states 
are present, we therefore find 


T,, = 


-(2D)(q ^ 0) + i = 0) 


(D5) 


Thus L~^T is the well-defined quantity in the thermody¬ 
namic limit, and only those edge-related terms that scale 
with L will contribute to it. 


Appendix E: Spin susceptibility of edge states 


with .=,2 (7) = 1 -bSy-b (7- 1)Clh and ^2,(7) = 1 -f (27- 
1)Clh = Ca;( 7 ). Thus oc where 

Xg®‘^ge)((7) is the Lindhard function associated with the 
edge states. Explicit calculation of (E3) yields 


xt^^\T,q) 


'='2 


(7) 


16t:A 

S^(7) 


327rA 


In 



, (E7a) 


-4fc|| \ 

4A2-g2 J ’ 


(E7b) 


where A is a large-wave-vector cut-off. Due to the spe¬ 
cial energy dispersion Eq. (El), q) is a constant 

(independent of kp and q). Clearly, for the hole doped 
case {fi < 0), the same result as in Eq. (E7) is obtained 
due to the assumed particle-hole symmetry. 


1. Spin susceptibility of edge states: Omitting 
finite size effects 


2. Spin susceptibility of edge states: Finite size 
effects included 


Following [5] (see also [44]), the dispersions for the edge 
states, using open boundary conditions for a confinement 
along the x direction and assuming D = 0, is given by 

e[^'> = sAk , (El) 

where k = ky. The associated spinor wave functions are 


2 

iQ)=Yl (^) V'of : (E2) 

where (/3Q^^(a:) = with cj)'^ = (l,±f) 

and Ai ^2 = ~^(1 i \/^ ~ 4|Cm|), since M < 0. To sim¬ 
plify matters, we have assumed particle-hole symmetry 
{D = 0) and the system length L —>• 00, i.e., the gap 
in the dispersions of the edge states is negligible. Note 
that the spinors in Eq. (E2) are independent of the wave 
vector component along the y direction. 

The spin susceptibility of the edge states is calculated 
by 


(edge) 


(7;'?; Q') 



(7; Q, q') 


where 

b') 


- nFlEk+l] 


(E3) 


[?7^*^(e)]^ ■ 4(7) • V^'’'Kb) 


x[4^4e0]^-4(7)-4^^(e')-(E4) 


With an averaging over the coordinates along the con¬ 
fined directions we have 

I de dg' B, b') = (J fj , (E5) 

JdedQ'w,^y\r,B,B') = ^(l J) , (E6) 


Following [44] , we now take into account the finite size 
of the system. As a result, edge states at the two sides 
that have the same spin can couple which results in a 
gapped spectrum of their energy dispersions.We assume 
the system size L to be large and use an approxima¬ 
tion of the wave functions in [44]. Taking into account 
the various degrees of freedom of both system sides, we 
modify the approach of [5], where the solutions for the 
wave functions read 

D _ ^ 

(E8) 


where C is a normalization constant and v = ± denotes 
the right- and left-mover. For non-zero wave vector fc, 
we consider the matrix 


"Hrl 


fAk A \ 
\A -Ak) 


(E9) 


in the basis of right- and left-mover, where A is the in¬ 
duced gap which is a function of the system parame¬ 
ters [44]. The energy dispersions of (E9) are 


= r\/(Afc)2 + A2 , (ElO) 


where r = ± labels the (spin-degenerate) conduction and 
valence bands, respectively. The associated eigenstates 
are 



(Ell) 


With this information, the spinors in Eq. (E2) are mod¬ 
ified to 

dir (^) = HI]«ir4 (x) (z) , (E12) 

1 — 1 V 



















12 


Thus, for the present case, the spin susceptibility of the 
edge states is 


(edge) 


(7;?; Q. q') 



2ti jj{k,k+q,T,T’) 


(7; Q, q') 


where 




(s') 


ihrj 


41-e: 


(s') 

k+qr' 


(E13) 


w(®"®') 

^jjik,k+q,T,T') 


(7; Q, Q')= 


[viliQ)1 ■ Sj{j) ■ villr'ie) 


x[vi+gr'ie'4 ■ Sjh) ■ vihe') ■ 


(E14) 

Averaging over the coordinates along the confined direc¬ 
tions, we obtain for the overlap factor of the Lindhard 
function 

Jdg = 

{ {akak+q + bkbk+qY (f = r') 

(E15) 

{(kkbk+q bkdk+q^ (f F ) 


The result for the overlap factors of the spin susceptibil¬ 
ities is 


J dg dg' 4 ltk\+q,r,r’lT. Q. q') 



a. Intrinsic contribution to the spin susceptibilities of the 
edge states in the limit g —>■ 0 


Calculating the intrinsic contribution to the spin sus¬ 
ceptibilities of the edge states in the long-wavelength 
limit (g —> 0), we obtain 


X 


(int,e) 

zz 


(7;0) 


(NLA)^ 


9sCUl) 

16fA 


1-ln 



(E19a) 


(int.e) 

Axx 


(7;0) 


gsSg(7) 

IGttA 


1 — In 



, (E19b) 


where A = A/A. We note that the intrinsic contribution 
to the Lindhard function vanishes in the limit g —>■ 0, 
which can be inferred from (E15). To compare this result 
with the one of the intrinsic bulk contribution, we divide 
(E19a) and (E19b) by the length L and let L go to infinity 
(see Sec. D). Thus, we obtain 


lim 

L—>oo 


L 


= 0 , 


(E20a) 


lim 

L—foo 


(int.e) 

Xxx 


(7;0) 


L 


9s-lh) 

327r|B| 


(1 - \/l - 4 |^m|) • 


(E20b) 


Therefore, the edge states give a contribution to the total 
susceptibility only for the in-plane component and its 
sign equals that of the bulk contribution. 


X 


and 


[(ufcUfc-l-g bkbk-\-q) 

+T{akbk+q -f bkak+q) Cz(n) NLAf 

[{bkbk+q - dkak+q) Czin) NLA 

^+T{akbk+q - bkdk+q) 2^(7)]^ 


(r = t') 

(E16) 
(t ^ r') 


dg dg' 4 liik+q,ry)(r, Q. q') = 

\dkbk+q bkdk-\-q (f — F ) 

+F(afeafc+, -I- bkbk+q)NLA\'^ 

[dkdk+q bkbk-\-q 
-T{akbk+q - bkak+q)NLAY 

In obtaining (E15)-(E17) we have used 


"2 


(7), 


X 


, . M (E17) 

(r ^ t') 


b. Electron doped contribution to the spin susceptibilities of 
the edge states in the limit g —>■ 0 


Next we include the contributions due to doping to the 
spin susceptibilities. The interband contributions read 
(for Akp A) 


X 


(inter,e) 
zz 


(7;0) 


-{NLAf 


9sCUl) 

IGttA 


1-ln 


(x)]' 

(E21a) 


X 


(inter,e) 

XX 


(7;0) 


98^1{i) 

IGttA 


1 — In 



(E21b) 


while the intraband contributions are 


x(^")(7;0) 


IGttA 


(E22a) 


f dxJJ « L , (E18) 

d u 

the functional L-dependence of A Ri F exp(A2L) [44], 
where FjEc, = 4|^m|/\/1 — 4|^m|, and we have defined 
N = (FC)-^. 


0) = -{NLAf , (E22b) 

which is consistent with the finding that 
{XxT'^^'’^\i',9)] are important [unimportant] for A —>■ 0, 
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whereas it is the other way around for the interband con¬ 
tributions. Considering 


lim 

L—¥C0 




(7;0) + x 


(inter,e) 
zz 


(7;0)] 


L 


0 , (E23a) 


lim 

L—¥oo 


L 


even when effects due to structural inversion asymmetry 
(SIA) are included. By taking into account the influence 
of a perpendicular electric field £z^ it has been shown in 
Ref. [22] that the BHZ Hamiltonian is supplemented by 
entries that mix the spin-up and spin-down components 
of the BHZ basis states. The leading contribution due to 
SIA is linear in the wave vector and given by 


f^(l-v'rV4M).(E23b) 

we find that this is the same contribution as the intrinsic 
contribution, Eq. (E20b), which has however the opposite 
sign. Thus, the sum of doped and intrinsic contributions 
of the edge states vanishes in the large L limit. As a 
consistency check, we verify for q = Q that the sum of 
Eqs. (E19), (E21) and (E22) yields Eqs. (E7a) and (E7b) 
(multiplied by gs) in the limit A —> 0. 


Appendix F: Effects of structural inversion 
asymmetry 

Here we demonstrate that the basic features of the in¬ 
trinsic spin susceptibilities given in Eq. (1) are robust 

I 


0^ 


/ 0 0 —iRok- 0\ 

0 0 0 0 

iRok^ 0 0 0 

V 0 0 0 0/ 


(El) 


where k± = kx i iky. The necessity to avoid dielectric 
breakdown provides an upper limit for the electric-field 
magnitude through the condition \e£z\d < 2|M|. Defin¬ 
ing the SIA-related dimensionless parameter = Rq/A, 
this condition translates into |^r| < 16.1 |^M|/(d[nm]) Ri 
0.12 for a typical heterostructure [22]. Thus is gen¬ 
erally a small parameter, and a perturbative treatment 
for SIA effects is appropriate. To lowest order in ^r, the 
intrinsic spin susceptibility in the limit g —> 0 is found as 


Xxx*n7;q 

xi“H 7 ;q 


cUi) 1 (.,.2 \ 

167rjHj I-b 4^M0('fM) V ^ 3[1-b 4 ^m0(Cm)]^ / 

CKx) 1 A , ,2 2(97^+ 2eM0(gM)[C.^(7) +187^]} 

167rjHj I-b 4^M0('fM) V ^ 3C2(7)[1-b 4 ^m0(Cm)]^ 


(E2a) 

(E2b) 


Thus the lowest-order SIA corrections to Xj™*^(7;q = 0) 
are quadratic in the small parameter ^r. This means that 
the result given in Eq. (F2) represents already an excel¬ 
lent approximation. Interestingly, Xxx^'^ (7; q = 0) turns 
out to not be modified by SIA contributions in the in¬ 
verted regime. We find that this remains true even when 
higher-order corrections in ^r are considered. In con¬ 
trast, Xz”*^(7; q = 0) in Eq. (E2b) has finite SIA correc¬ 


tions in the inverted regime given by Sy^^R/C^ (7). In our 
case where 7 = —2.22 this amounts to a relative change 
that is about 1%. Also in the normal regime, SIA contri¬ 
butions to (7; q = 0) are at most of relative magni¬ 
tude 1%. Thus we conclude that the spin susceptibilities 
in Eq. (1) generally receive only very small corrections 
when SIA terms are included in the BHZ Hamiltonian. 
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